Far-from-equilibrium Ostwald ripening in electrostatically driven granular powders 



M.V. Sapozhnikov 1 ' 2 , A. Peleg 3 , B. Meerson 4 , LS. Aranson 1 , and K.L. Kohlstedt 5 
1 Argonne National Laboratory, 9700 S. Cass Avenue, Argonne, IL 60439 
2 Institute for Physics of Micro structures, Russian Academy of Sciences, GSP-105, Nizhny Novgorod 603600, Russia 
3 Theoretical Division, Los Alamos National Laboratory, Los Alamos, NM 87545 
A Racah Institute of Physics, Hebrew University of Jerusalem, Jerusalem 91904, Israel and 
5 Department of Physics and Astronomy, University of Kansas, Lawrence, KS 66045 

We report the first experimental study of cluster size distributions in electrostatically driven 
granular submonolayers. The cluster size distribution in this far-from-equilibrium process exhibits 
dynamic scaling behavior characteristic of the (nearly equilibrium) Ostwald ripening, controlled by 
the attachment and detachment of the "gas" particles. The scaled size distribution, however, is 
different from the classical Wagner distribution obtained in the limit of a vanishingly small area 
fraction of the clusters. A much better agreement is found with the theory of Conti et al. [Phys. 
Rev. E 65, 046117 (2002)] which accounts for the cluster merger. 

PACS numbers: 45.70.Mg, 45.70.Qj 



Dynamics of big assemblies of macroscopic particles 
are poorly understood, especially when inter-particle in- 
teractions are strongly dissipative [J, Q . Additional com- 
plications arise when the grain size goes below 0.1 mm, 
and non-trivial contact interactions come into play. As 
small particles acquire an electric charge, the dynamics 
are governed by the interplay between long-range elec- 
trostatic and short-range contact forces. On the other 
hand, an efficient electrostatic excitation of granular me- 
dia becomes possible here, and offers new opportunities 
for investigation and control. 

Recently, a host of fascinating pattern-formation phe- 
nomena and far-from-equilibrium phase separation pro- 
cesses have been observed in electrostatically driven gran- 
ular powders, in air or vacuum 0, 0, 0, an d in poorly 
conducting liquids 7]. One of them is cluster coarsen- 
ing which is remarkably similar to the Ostwald ripen- 
ing H E|, despite the fact that the system is very far 
from equilibrium. This work deals with experiment and 
theory of this process. We work with submonolayers of 
metallic powders, i.e., the layer of particles only partially 
covers the bottom of the experimental cell, and it is ap- 
proximately one-particle high. Conducting particles (we 
used 40 fim spherical copper particles, the total number 
of particles was about 10 7 ) are placed between the plates 
of a large plane capacitor, levelled horizontally. We used 
27 x 27 cm transparent capacitor plates (glass coated by 
indium doped tin oxide); the plate spacing was 1.5 mm. 
The DC electric field E in the capacitor was varied in 
the range — 7 kV/cm. The experiments were performed 
in the atmosphere of dry nitrogen to reduce adhesion of 
particles on the plate due to humidity of air. Visualiza- 
tion was achieved with a CCD camera mounted above 
the cell. 

When particles are in contact with the capacitor plate, 
they acquire an electric charge. If the electric field in the 
cell E exceeds a critical value, the resulting upward elec- 
tric force overcomes the particle weight and pushes the 
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FIG. 1: A schematic of experimental setup. 



particles upward. When a particle hits the upper plate, 
it recharges and falls back. Remarkably, the collective 
behavior of this far-from-equilibrium system closely re- 
sembles that of phase separating systems close to equi- 
librium 0, The role of temperature, however, 
is played by the electric field E. The particles remain 
immobile on the bottom plate of the capacitor, and form 
a precipitate phase, at E < E\. If the field E is larger 
than a second threshold value, E2 > E±, the particles 
form a gas-like phase. This second field E2 is 50%-70% 
larger than E\. Upon "quenching" the field into the co- 
existence region E\ < E < E2 one observes nucleation of 
precipitate, and small densely packed clusters form and 
start to grow on the bottom plate. Then the clusters 
undergo coarsening: larger clusters grow at the expense 
of smaller ones (which shrink and disappear). This far- 
from-equilibrium coarsening process closely resembles the 
Ostwald ripening with attachment- detachment (in a dif- 
ferent terminology interface- controlled) kinetics |9(. 

The basic physics behind phase separation and coars- 
ening in this system was investigated previously 0, IS • 
Two main physical processes are at work here. On the 
one hand, there is electrostatic screening: a decrease in 
the vertical electric force F, exerted on a grain in contact 
with the bottom plate, caused by the presence of other 
grains. On the other hand, collisions between the grains 
and the bottom plate, and between the grains themselves, 
are strongly inelastic, so the granular temperature is of 
no importance. The reader is referred to Refs. 0, for 
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detail, including quantitative estimates of the respective 
forces etc. 

At the coarsening stage the clusters compete for the 
material, and this competition is mediated by the gas 
phase. A preliminary scaling analysis of experimental 
data on the cluster coarsening dynamics Q and special 
measurements of the shrinkage rate of single clusters Q 
strongly indicated that the cluster growth/coarsening is 
limited by the attachment-detachment rates of the "gas" 
particles, rather than by the transport rate in the "gas" 
phase. Correspondingly, a phenomenological continuum 
model of this phase separation and coarsening was formu- 
lated in Ref. p in terms of a Ginzburg-Landau equation 
subject to conservation of the total number of grains. 
Similar globally-conserved Ginzburg-Landau equations 
have been studied in different contexts [l(J, |l2| . In 
the regime of well-developed (densely-packed) clusters 
the continuum model yields "sharp-interface" equations 
which govern the dynamics of the inter-phase boundary 
and are very useful in the analysis of the late-time stage 
of the phase separation 

This work focuses on the statistical dynamics of this 
far-from-equilibrium Ostwald ripening process. One ad- 
vantage of our new experimental apparatus is a very 
large aspect ratio, allowing us, for the first time, to pro- 
duce a large number of clusters and accumulate sufficient 
statistics. By following the dynamic scaling properties 
of the cluster size distribution we establish with confi- 
dence the attachment-detachment character of kinetics. 
Surprisingly, the scaled size distribution of the clusters 
turns out to be quite different from the classical Wagner 
distribution A much better agreement with exper- 
iment is found when we employ a more advanced the- 
ory by Conti et al. which, in addition to the clus- 
ter growth/shrinking processes, accounts for the cluster 
merger. The theory has no adjustable parameters and is 
very different, both in its assumptions and predictions, 
from a recent phenomenological exchange-driven growth 
model suggested by Ben-Nairn and Krapivsky 13]. 

Figure 2 shows a typical sequence of snapshots of the 
system in the coexistence region E\ < E < E2. The 
Ostwald ripening process is clearly seen. In addition, 
distinct merger events are observed. Figure 3a shows 
the average cluster area (A)(t) versus time. This de- 
pendence is very close to linear. Introducing the cluster 
radius R = (A/ir) 1 / 2 , we find that the mean-square ra- 
dius R*(t) = (R 2 ) 1 / 2 grows with time like t 1 ! 2 . This is a 
defining feature of the attachment-detachment-controlled 
Ostwald ripening 0|. The plot of the inverse number 
of clusters in the cell N~ x versus time (Fig. 3b) ap- 
proaches a straight line which implies a scaling behavior 
N(t) ~ t~ x . An additional support for the scaling is pro- 
vided by the approximate constancy of the total area of 
the clusters at late time, N(t)R*(t) 2 ~ const, see Fig. 
4. The dynamic range of Figs. 3 and 4 is limited by 
5.40 x 10 4 , since for later times the number of clusters is 
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FIG. 2: Snapshots of coarsening of granular clusters in elec- 
trostatic cell at times (in seconds) t — (upper left), 1 x 10 4 
(upper right), 2 x 10 4 (lower left) and 5 x 10 4 sec (lower right). 
The applied DC electric field E = 2.33 kV/cm. 
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Fit: a+bt c 




a 23.6 ±1.7 




b 0.60 ±0.05 




c 1.00 ±0.014 
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Fit: a+bt c 




- a 0.001 65±9.0E-5 




: b 3E-5 ±2.3E-6 




c 1.042 ±0.014 






(b) 



10 



100 



FIG. 3: The average cluster area (A) (a) and the inverse 
number of clusters in the cell 1/N (b) versus time in the ex- 
periment shown in Fig. 2. The unit of time in this figure 
corresponds to 100 sec. The solid curves (see also the insets) 
show data fits of the type a + bt c . The straight line in figure 
b, shown for reference, has slope 1. Notice that in the exper- 
iments we waited for about two minutes (until small clusters 
nucleated at the bottom plate) before starting the measure- 
ments. As a result, the (A)(t) curve does not pass through 
the origin. 



too small. 

A more detailed characterization of the far-from- 
equilibrium Ostwald ripening is provided by the dynam- 
ics of the probability distribution function (PDF) of the 
cluster sizes f(R,t). The PDF is defined as usual: 
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f(R,t) dR = N(t)/L 2 



(1) 



where L 2 is the cell area. Dynamic scaling behavior of 
f(R,t) is expected at late times, when the total area of 
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FIG. 4: The area fraction of the clusters (the total area of the 
clusters divided by the area of the cell) versus time. The unit 
of time corresponds to 100 sec. The dynamic scaling regime 
is observed at late times, when the area fraction approaches 
a constant value. 



the clusters approaches a constant: 

POO 

7T / f(R,t)R 2 dR = e. 
Jo 



(2) 



Here e is the area fraction of the clusters: the ratio of 
the total area of all clusters to the cell area L 2 . In view 
of Eq. J2J), the expected scaling behavior is 



f(R,t) = R: 3 (t)F[R/R*(t)] 



(3) 



Using the Ansatz (PJ), we can express the number of clus- 
ters JD an d the mean squared radius R%(t) through the 
moments of the scaled PDF F(^): 



N(t) 
L 2 



m 
Rl{t) 



(4) 



Kit) = (R 2 ) 



_ f™f(R,t)R 2 dR _ m 2 



J™f{R,t)dR 



m 



Rl(t), (5) 



where rrik = J °° F(£) £ k dt; , k = 0,1, Equation © 

yields a consistency condition rri2 = mo. 

The imaging processing and tracking of the cluster 
growth and coarsening in our experiment was performed 
using the MATLAB 6 toolbox. The scaled PDF, obtained 
by averaging over 4 different experiments, performed un- 
der the same conditions, is shown in Fig. 5. The hori- 
zontal axis shows the cluster radius divided by the 
vertical axis shows the experimental PDFs multiplied by 
Rl(t). The scaled PDF is the result of collapse of 340 
original, unsealed PDFs from the 4 experiments, mea- 
sured in the time range 2.01 x 10 4 sec < t < 5.40 x 10 4 
sec with equal 100 sec time intervals. We checked that 
the late-time values of the cluster area fractions in these 
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FIG. 5: Scaled probability distribution function (PDF) of 
cluster sizes where f = R/R*(t) [see Eq. ©]. The 

squares show the experimental results, obtained from 340 
snapshots from the interval 2.01 x 10 4 sec < t < 5.40 x 10 4 
sec and averaged over 4 experiments. The error bars are not 
shown when they are smaller than the size of the square. The 
dotted line is the Wagner distribution Fw(£) [see Eq. @], for 
the same area fraction e — 0.092 as in the experiment. The 
solid line shows the scaled PDF, obtained from the theory of 
the attachment-detachment-controlled Ostwald ripening with 
coalescence [3| f° r the same area fraction. 



4 experiments were close to each other: 9.1%, 9.3% (in 
two experiments) and 9.4%. 

What is the theoretical prediction for the scaled PDF 
F(£)? Let us first proceed by neglecting the cluster 
merger. In this case the PDF dynamics is described by 
a continuity equation in the space of cluster sizes: 



dtf + d R (Rf) = 0. 

For the attachment-detachment-controlled kinetics, 
cluster growth/shrinking law is 



R = D 



Rc(t) 



(6) 
the 

(7) 



where R c (t) is the time-dependent critical radius, while 
the effective diffusivity D depends on the electric field 
E and on the size and weight of the grains 5]. The 
growth law J7J) is obtained under assumption that the 
particle transport in the gas phase is very fast, so that 
the concentration of the gas phase in the cell is approx- 
imately uniform in space, varying only in time. The in- 
verse critical radius 1/R c (t) in Eq. J7J) is proportional to 
the (time-dependent) difference between the concentra- 
tion of the gas phase and a constant value of this concen- 
tration, for which a planar interface of the cluster phase 
is at rest 5j. Equation Q describes shrinking and dis- 
appearance of small clusters, for which R < R c . The 
shrinking occurs in a finite time. At R <C R c the area 
of the shrinking cluster goes down linearly as the time 
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increases. This property was checked in experiment, and 
was found to hold until the shrinking clusters become 
small compared to the gap between the capacitor plates, 
and three-dimensional effects come into play @]. 

Equations (0 and Q, together with the cluster 
area conservation 0, make a closed set. Mathemati- 
cally, this model is identical to the Wagner's model of 
the attachment-detachment-controlled Ostwald ripening 
0. (The attachment-detachment-controlled systems are 
quite different from the diffusion-controlled systems, for 
which a theory of cluster coarsening, valid for vanishingly 
small area fractions of the clusters, was developed earlier 
by Lifshitz and Slyozov An immediate consequence 
of Eqs. 0, ® and (0) is the equality 



R c (t) = (R)(t) 



f~f(R,t)RdR 
f °°f(R,t)dR 



For the scaling solutions © we obtain R c (t) = (R)(t) = 
(mi/mo)R*(t). Therefore, at large times (R)(t) = 
const (Dt) 1 / 2 and, in view of Eq. N(t)/L 2 = 

const (Dt) -1 , where the prefactors depend on e. Em- 
ploying the Ansatz (0 and solving the resulting ordinary 
differential equation, one arrives at a family of solutions 
for all of which have a compact support < £ < £ m . 
The problem of selection of the correct scaled PDF out 
of this family of solutions has been extensively studied 
[ill Il4l H5I ll^ . A special role here is played by the 
Wagner distribution, 



Fw(0 




exp 1 



?) ifO<£<£ m , 

if £ > fm, 



(8) 



-1/2 



for which £ m = [l + 2e 2 Ei(-2)] ' ~ 1.8989 is the 

largest. Also, C = tt' 1 [l/(2e 2 ) + £i(-2)] _1 ~ 16.961, 
and Ei (...) is the exponential integral function [Tjj. The 
Wagner distribution for e = 0.092 (the same area fraction 
as in experiment) is shown in Fig. 5. One can see a large 
discrepancy between the Wagner distribution and the one 
observed in experiment. 

The Wagner's model is inaccurate because it disregards 
the cluster merger. Therefore , we employed the more ad- 
vanced theory of Conti et al p2 which accounts, within 
the framework of the Ostwald ripening, for the binary co- 
alescence of clusters (growing clusters merge upon touch- 
ing each other). Now Eq. © gives way to the following 
kinetic equation: 



1 r 

2 Jo 



d t f + d R (Rf) 



6(R - #1) + 5(R-R 2 )-5[R- J Rj 



{2M(R U R 2 ) 



R\ 



f(R 1 ,t)f(R 2 ,t) dR 1 dR 2 },(9) 

where M(R U R 2 ) = tt (R ± + R 2 ) (Ri + R 2 ) 0(Ri + R 2 ), 
6(. . . ) is the step function, and R\ and R 2 are gov- 
erned by the same growth law J7J). We refer the reader 



to Ref. f° r a derivation of Eq. ©. Remarkably, 
Eqs. (0), and © admit the same scaling Ansatz 
(JHJ). The scaled PDF is described by a nonlinear integro- 
differential equation [which follows from Eq. Q], subject 
to m 2 = mo = e/n. The integro-differential equation can 
be reduced to an integral equation and solved iteratively 
[l^ . A crucial role in the iteration procedure is played 
by the (non-normalizable) solutions of the Wagner prob- 
lem which have an infinite support < £ < 00. The 
solution of the problem is unique (and normalizable), so 
the account of merger resolves the selection problem in- 
trinsic to the Wagner's formulation. Figure 5 shows the 
scaled PDF, obtained by the iteration procedure [l2| for 
the same area fraction e = 0.092 as in the experiment. 
One can see that the agreement is much better than for 
the Wagner distribution. Importantly, the only parame- 
ter which enters the theory of the scaled PDF is the area 
fraction e. Therefore, for a given area fraction, a com- 
parison of the theory and experiment does not involve 
any adjustable parameters. As can be seen from Fig. 5, 
the theory works very well in describing the position of 
the maximum of the scaled PDF and the shape of the 
PDF to the right of the maximum, including the tail. 
On the other hand, the theory overestimates the number 
of clusters at very small scaled radii and underestimates 
it in the region of the maximum of the PDF. This rela- 
tively small but systematic disagreement may result from 
the assumption of an infinitely fast transport in the gas 
phase. If this assumption is relaxed, inter-cluster cor- 
relations will appear which should reduce the effective 
merger rate. In our experiment, however, this reduction 
is mitigated by the fact that, as time elapses, more ma- 
terial gets concentrated in the middle of the cell (see Fig. 
2). The future work should attempt to account for these 
two competing effects. 

In summary, we investigated the statistical dynam- 
ics of clusters in electrostatically-driven granular pow- 
ders. We found that a recent advanced theory of the 
attachment-detachment-controlled Ostwald ripening [l2[ 
yields a good quantitative description of this far-from- 
equilibrium system. It is tempting to also apply this kind 
of analysis to mechanically vibrated granular monolay- 
ers which exhibit strik ingly s imilar non-equilibrium phase 
separation properties 0, LU| • 
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